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ABSTRACT 



Aims. We analyze the ratios of the broad hydrogen Balmer emission lines (from Ha to He) in the context of estimating the physical 
conditions in the broad line region (BLR) of active galactic nuclei (AGN). 

Methods. The Balmer emission lines are obtained in three ways: i) using photoionization models obtained by a spectral synthesis code 
CLOUDY; ii) calculated using the recombination theory f or hydrogenic ions; iii) measured from the sample of observed spectra taken 
from the Sloan Digital Sky Survey database. We investigate the Balmer line ratios in the frame of the so called Boltzmann-plot (BP), 
analyzing physical conditions of the emitting plasma for which we could use the BP method. The BP considers the ratio of Balmer 
lines normalized to the atomic data of the corresponding line transition, and is in that way different from the Balmer decrement. 
Results. We found that for a certain range of thermodynamic parameters, there are objects that follow the BP. These AGN may have 
a BLR consisting of mostly high density plasma. 

Key words, galaxies: active - (galaxies:) quasars: emission lines - line: formation - plasmas 



1. Introduction 

The investigations of the broad line region (BLR) structure (size, 
geometry, physics, etc.) of active galactic nuclei (AGN) are 
nowadays still important. One of the main implications is the 
accurate estimate of the mass of the super-massive black hole 
(SMBH) in the center of an AGN, since it can be derived from 
the dyna mics of BLR gas, gravitationally bound to the SMBH 
(see e.g. iMcGill et all I20081 iBentz et all l2010h . After several 
decades of research, our view of the physics and structure of 
the BLR is improving, b ut we sti l l lack a full understanding 
of its properti es (see e.g. iBaldwinl Il997b ISulentic et all l2000t 
lGaskelll2009l for a review). It is obvious from the broad emis- 
sion line (BEL) fluxes and profiles that the BLR physics and 
kinematics are more complicated than in the narr ow line region 
(NLR) or oth e r nebular environments (see e.g. | Kroli5 |1999t 
ISulentic et ail l2000l lOsterbrock & Ferlandi l2006l and refer- 
ences therein), and they are even closer t o the conditions in stel- 
lar atm ospheres than in gaseous nebulae (lOsterbrock & Ferlandi 
I2006I) . 

Most of the BLR gas is photoionized. This was first inferred 
by the approxim ately constant value of the H/3 equivalent width 
in several AGN (lOsterbrockl Il978l) and then confirmed by re- 
verberation studies showing a direct res ponse of the BEL fluxes 
to the continuum flux changes (see e.g. iPetersonl 1 1 9941) . So far 
there have been numerous attempts to describe the BLR physics 
and to explain the emis sion-line spectrum by means of photoion- 
ization models (see e.g . lNetzetl 19751 Davidson & NetzeH l 19791 
Kwan & Krolilll98lll^llin-Souffrin. Dumont & Tullvl 19821 
Rees. Netzer & Ferlandi Il989l iBaldwin et all Il995l Il996 : 
Dumont, Collin-Souffrin & Nazaroval I1998T | Ferland et al. 
19981 iKrolikl 1 19991 iMarziani et ail l200lE iFerland . 



20031 iKorista & Goadl l 200l iLeighlv & Casebeeil [20071: 
Snedden & Gaskelll l2007l IMarziani et all l2010h . Today there 



are roughly half a dozen major c odes, that were continu ously 
developed since the 1970s (see e.g. lFerland & Savi nl l200Tl for a 
review). Although such numerical simulations largely improved 
in the course of years, the investigations of the BLR still contai n 
some important unresolved problems (see e.g. iFerlandl 119991) . 
One of the most interesting concerns the observed small ratio 
among the Lya and H/3 lines (usually in the range of 5 - 15), 
while mode l s pred ict much higher values (30 - 50) (see e.g. 
iNetzer et all 1 19951 and references therein). Another important 
finding is that there are indications for a physical complexity of 
the BLR, which could result from the combination of distinct 
components, as suggested, for instan ce, by the origin of "high" 
and "low ionization l ines" (see e.g . ICollin-Souffrin & Lasotal 
119881 IBaldwinl 1 19971 ISulentic et all 120001) . but also according 
to the complex profiles of single lines, that cann ot be explained 
in terms of a simple single-region model (see e.g . Popovic et all 
12001 IPopovicl l2006bl: iBon et all 120091 [Gaskell l2009l )~ 

Different BELs are originated at various distances from the 
central continuum s ource and under a w ide range of physical 
conditions (see e.g. ISulentic et all l2000t) . The BLR is of high 
density, thus the collisional excitation, self-absorption, dust ob- 
scuration and complicated coupling of line and continuum ra- 
diation transfer should be taken into account in the calcula- 
tion of the resulting spectrum (see e.g. iBaldwinl Il997l) . In such 
circumstances, the general spectroscopic techniques, commonly 
adopted to estimate physical conditions in nebular environments, 
cannot be applied. Some clues, such as the observed Fe II emis- 
sion, suggest that T < 35,000 K, since at higher tempera- 
tures it would be effectively suppressed by collisional ioniza- 
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tion to Fe III (lOsterbrock & Ferlandl 12006). However it is very 
likely th at the Fe II emission only affects a fract i on of the BLR 
(see e.g. iMarziani et all 120011: lllic et all 120091 : IPopovic et all 
I2009L iKova cevic et alll2010l and references therein). The elec- 
tron density is estimated to be in the range 10 8 cm~ 3 < n e < 
10 14 crrT 3 , in order to suppress the emission of broad forbidden 
lines, while still allowing for the presence of p ermitted and semi- 
forbidden ones (lOsterbrock & Ferlandl 120061) . For example, for 
the upper limit the CIII] /U909 emission line implies that the 
densities cannot be higher than ~ 10 n ~ 13 cirT 3 , though strong 
UV lines like Fe II or Al III A 1860 suggest somewhat higher den- 
sity, at least in the low ionization reg ion (lBaldwinetallll996t 
lLaor et allfT997tlNegrete et all 1201 lh . 

There are not many proposed methods in the literature that 
include obser vations of the BELs to determine t he BLR physi- 
cal pro perties (IMarziani et all 1200 it lLaorl 120061) . IMarziani etaT] 
( 12001b found, using the CLOUDY photoionization computation, 
that the ratio of Si III] ^1892 and C III] A 1909 UV emission 
lines is a good density diagnostic in the density range 9.5 $ 
logn e S> 12. And since that ratio is correlated with the width 
of the broad H/3 line, they gave a relationship for the estimates 
of the electron density in the BLR using either the ratio of these 
semi-fo rbidden lines or the wi dth of the broad component of the 
H/3 line (IMarziani et al 1 120011) . The same authors extended their 
analysis and analyzed the physical conditions using other line 
ratios (e.g. Al III ^1860/Si III]^1892, Si III]^1892/C IV^1549, 
Mg IU2800/Lyff, etc.) trying to exploit the CLOUDY simula- 
tions to deduce c onstraints on ionizatio n para meter, densit y and 
column density ( IMarziani et al 1 l2010h . Also, lLao n (120061) pro- 
posed a method that considers the electron scattering influence 
on the line profiles and determines the physical parameters of 
the BLR in the case of low luminosity AGN with the emission- 
line profiles having exponential wings (e.g. the case of galaxy 
NGC 4395). The method assumes that the exponential wings are 
produced by the optically thin, isotropic, thermal electron scat- 
tering. In that case, by fitting the line wings with an electron- 
scattering model, one ca n estimate th e electron density and op- 
tical depth of the region £aorll2006l) . Both described methods 
are observationally constrained for either needing the UV obser- 
vations or detecting BEL with strong exponential line wings to 
extract a direct estimate. 

On the other hand, IPopovic et al] (120021) : iPopovid (120031 
l2006al) suggested that the Boltzmann-plot (BP) method (see §2 
for more details), already well known to laboratory diagnos- 
tics of high density plasma dGrieml Il997l) . m ight be exploited 
to probe the BLR of some AGN (see a lso lllic et all l2006t 
lLa Mura et"ail 120071: IPopovic et all 120081) . By measuring the 
flux of emission lines belonging to a specific transition series, 
such as the hyd rogen Balmer line series in the optical domain 
dPopovid 120031 and using the corresponding atomic parame- 
ters, the technique provides an estimate of the plasma tempera- 
ture. For example, in case of NGC 5548 the BLR temperature 
was estimated usin g the BP applied to th e Balmer lines observed 
from 1996 to 2004 iPopovic et alll2008l) . A high correlation be- 
tween the variation o f the optical AGN co ntinuum and the BP 
temperature is found (IPopovic et al.L 12008b . 

The main motivation of this paper is to investigate under 
which particular circumstances the Balmer line ratios and BP 
method may, or may not, be used to explore the physical con- 
ditions of the BLR, given that photoionization is expected to 
control the plasma physical properties. To clarify this, we com- 
pare the hydrogen Balmer lines (from Ha up to He), that we 
obtained in three different ways: i) from a grid of several numer- 
ical models, computed using the CLOUDY spectral synthesis 



code dFerland et all Il998l) and analyzing the properties of the 
resulting spectra; ii) by considering the calculated emission line 
inte nsities predicted in the fram ework of the recombination the- 
ory (IStorev & Hummerl Il995al) . and iii) analyzing a set of ob- 
served sp ectra taken from the S loan Digital Sky Survey (SDSS) 
database ( lLa Mura et alll2007l) . 

The paper is organized as follows: in §2 we present the ap- 
plication of the Balmer line ratios through the Balmer decrement 
and BP method, in §3 we describe the grids of models generated 
with the CLOUDY code and report the results of the models, and 
we describe and analyze the theoretical recombination emission 
line, in §4 we report on the observed SDSS spectra, in §5 we 
provide discussion and finally, in §6 our conclusions are given. 



2. The ratios of the hydrogen Balmer lines - 
indications of the BLR physical properties 

The hydrogen Balmer lines are the brightest recombination lines 
in the optical spectra of AGN. They are usually showing com- 
plex line profiles, composed from the narrow and broad compo- 
nent. Here we discuss only the broad component, coming from 
the BLR. 



2. 1 . The Balmer decrement 

The hydrogen Balmer decrement - Ha /Hp - ratio 
used to probe the physic s of the line-emitting 



used to probe me p 
(lKwan&KrolikL[T97^ 
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plasma 



1981; Mathews, Blumenthal & Grand 



19801: ICanfield&Puetten, ll98lUDong et all 120081: iJin et al 



2012). It is most frequently used to determine the amount of dust 



extinction for th e low-density gas, such as the narrow line region 
(NLR) in AGN lOsterbrock & Ferlandl (120061) . where a value of 
3.1 for Ha/Hp ratio is generally adopted. This is slightly higher 
than the standard case B recombination value as the Ha emission 
is en hanced due to collisional effects and harder ionizing cont in- 
uum lGaskell & Ferlandl (ll984l) : IOsterbrock & Ferlandl (I2OO6I) . 

The densities in the BLR are far higher, as discussed above, 
so that the Ha/Hp ratio might be affected by high-d ensity effects 
or wa v elength-dependen t extinction by dust (see e.g. lOsterbrockt 
119841: iGoodrichl 11995b . The observed broad-line Ha/H/3 ra- 
tio of most broad-line AGN is usually la rger (steepe r ) than 
the Case B recombination value (see e.g. lOsterbrockl Il977l 
lRafanelh1 . ll985l;rDong et aUl2005l). sometimes even as steep as 
10 or h igher (|Osterbrock| 1198 It ICrenshaw. Peterson & Wagnerl 
Il988l) . iDong et alJ (l2008h obtained the averaged value of the 
broad-line Balmer decrement of 3.06, for a large, homogeneous 
sample of 446 low-redshift (z<0.35) blue type 1 AGN with the 
minimal dust extinction effects. They argue that the broad-line 
HafHfi ratio can be used as dust extinction estimator, even fo r 
the BLR, especially for the radio-quiet AGN dDong et alll2008l) . 

On the other hand, in several monitoring campaigns, the 
Balmer decrement is found to vary, us ually anti-correlated with 
the continuum flux, in a sin gle object (IShapovalova et all 120041 
l2010t IPopovic et all 1201 ll) . but often exhibiting complicated 
behavior depending on the state of the activity in the AGN. 
In case of NGC 4151, it varied from 2 to 8 during 11 years 
(IShapovalova et all 120081) . usually decreasing when the contin- 
uum increases, but staying constant in case of high continuum 
flux. A similar effect was detected in the case of 3C 390.3: for the 
low activity period the Balmer decrement is anti-correlating with 
the continuum, while in the outbursting phase with higher values 
of the c ontinuum flux, it stay s almost constant with the value 
of ~4.5 (IPopovic et all 1201 ll) . Thus the changes in the Balmer 
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decrement might indicate that sometimes the line production is 
not fully governed by the input ionizing continuum coming from 
the accretion disk, i.e., the main cause of the Balmer decrement 
variations is not related to the active nucleus and the shock ini - 
tiated emission is probably dominant (IShapovalova et al 
Or in addition to the physical conditions across the emitting re- 
gion or input ionizing flux, the size of the em itting regions could 
play an important role dPopovic et al.Ll20lTl) . Concerning the in- 
trinsic host galaxy reddening in these galaxies, we might expect 
that during the considered time interval, the intrinsic reddening 
can be neglected, because it should not vary too much for such a 
relatively short period (~10 years). Of course, the effects of dust 
within the BLR might still be present. 

As we noted above, the Balmer decrement (only the Ha/HB 
ratio) may indicate some physical processes in the BLR, but it 
is better to take into account the ratios of more lines from the 
Balmer series. One way to consider ratios of several lines is 
to use the Boltzmann-plot. We give a short outline of the BP 
method in the next Section. 



2.2. The Boltzmann-plot 

If a plasma extending over a region of length I emits along the 
line of sight, assuming that the temperature and emitter density 
do not vary too much, the flux F u i of a transitio n from an up- 
per level u to the lower I can be calculated as (lGriemill997t 
lKonievid[l^99llPopovid.[2003Ll2006al:[Popovic et alll2008l) : 



he N 
F u i = — guKi— exp(-E u /kT)£, 

A Z 



(1) 



where A is the transition wavelength, g u is the statistical weight 
of the upper level, A„; is transition probability, No is the aver- 
aged total number density of radiating species which effectively 
contribute to the line flux without being absorbed, Z is the par- 
tition function, E„ is the energy of the upper level, T is the av- 
eraged excitation temperature, and h, c, and k are the Planck 
constant, the speed of light, and the Boltzmann constant, respec- 
tively. Eq. (1) additionally assumes that the occupation number 
of the upper level in the transition follows the Saha-Boltzmann 
distribution, the lines are optically thin, and the line series orig- 
inates in the same emitting region (for more detailed derivation 
and discussion on th ese assumptions see iPopovicl 120031 l2006bt 
iPopovicet aill2008l) . 

Introducing a "normalized" line-flux F„, with respect to the 
atomic constants which characterize the transition: 



F n = 



Fui ■ A 



g U A u l ' 

we see from Eq. (1) that it is possible to write: 
log 10 (F n )=B-AE„, 



(2) 



(3) 



where, considering lines that belong to a specific transition se- 
ries, B and A represent the Boltzmann-plot parameters, with 
B being a characteristic constant of the transition series and 
A = log 10 (e) /kT exc x 5040/r exc is called the temperature pa- 
rameter. 

Therefore, taking into account a specific line series, such as 
the Balmer lines, if the populations of the upper energy states 
(n > 3 13 follow a Saha-Boltzmann distribution, Eq. (3) pro- 
vides an estimate of the average excitation temperature in the 
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Fig. 1. Examples of AGN where the BP method works (upper 
panel) and where it does not give valid results (bo ttom panel) 
applied on the data taken from lLa Mura et al ]d2007|).Each panel 
contains the observed spectrum and the BP applied on the line 
fluxes of the Balmer series (only the broad component of the 
Balmer lines was used). For the lower plots, which give the ob- 
served spectra, the observed wavelength scale (x-axis) is con- 
verted into energy units (eV), according to the energy above the 
ground of the upper transition levels, to match the BP in the up- 
per plots. 



region where the lines are formed. In the case of a low den- 
sity environment, where the nar row lines come from, the above 
assumptions do not apply (see iPopovid, 120031) . but for higher 
densities in the BLR plasma, the approximation might be rea- 
sonable. The Balmer line ratios are less sensit ive to temperature 
changes at higher temperatures (iGriem 1 H997I) . For further anal- 
ysis we adopt the maximal value of the BP temperature to be 
30,000 K, as the BLR temp erature should be lower than that 
dOsterbrock & Ferlandi 120061) . 

The Boltzmann-plot should not be confused with the Balmer 
decrement, as it considers line-ratio normalized using the atomic 
parameters of the specific line transitions. However, all discus- 
sion presented for the Balmer decrement in the previous section, 
applies for the BP since the line ratios could be highly influ- 
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Fig. 2. Different examples of the BP applied to the Balmer lines 
obtained with the CLOUDY models. The intensities F n are cal- 
culated using the Balmer lines normalized to the H/? flux. In the 
upper right corner, the BP temperature Tbp and the error of the 
BP fit / (in percentage) are given. In the bottom left corner the 
column density A^h, the hydrogen density «h and the input ion- 
izing flux Oh are given. 



enced by the dust effects. On the other hand, the BP method 
has the clear advantage to require only the measured Balmer 
line fluxes to estimate the excitation temperature. But one has to 
consider some possible drawbacks, concerning the use of emis- 
sion lines to infer the BLR physical properties, which in general 
appear in all methods that use BEL parameters for plasma diag- 
nostics. Since the BEL profiles are complex, the estimate of their 
fluxes should deal with the possibility of a multiple-component 
structure in the BLR, as stated in previous section. Furthermore, 
there are some indications suggesting that the Balmer lines do 
not necessarily have to arise strictly from the same region, as it 
is point ed out by some differences in the profiles of H a and HJ3 
(see e.g lShapovalova et all 120081 lPopovi6 et alll201 lh . Finally, 
different mechanisms may contribute to the formation of the 
emission lines and, though photoionization is the main heating 
source, other processes might be relevant. The combination of 
these effects implies that the actual properties of the BLR may 
or may not fit in the fundamental assumptions of the BP method, 
as it is illustrated in the examples plotted in Fig.Q] 



3. Simulated and theoretical Balmer line ratios 

3. 1 . Simulation of photoionized region and modeled Balmer 
lines 

Numerical simulations make it possible to understand complex 
physical environments starting from fundamental principles of 
physics. This is particularly applicable in the case of the BLR 
where due to different physical conditions in the plasma (espe- 
cially high density), many physical processes should be taken 
into account. To model such a complex physical s cenario, our 
choice was to use CLOUD YFI dFerland et all fl998h . Collisional 
effects, including excitation and de-excitation, continuum fluo- 
rescence, recombination, etc. , are al l accounted for in this code 
dFerland et all [l99i IFerlandl l2006h . 

In order to investigate the Balmer lines, we generate grids of 
photoioniz ation models of the BLR using the v ersion C08.00 of 
CLOUDY dFerland et all Il998l IFerlandl l2006h . This version is 
roughly a factor of 2 slower than C07.02, due to a major expan- 
sion in the physics of the H-like iso-sequence. The 1-levels of 
the H-like sequence are now fully resolved, making these atoms 
formally correct, thus the predicted H I emission line spectrum 
is changed with respect to the C07.02. Input parameters for the 
simul ations are chosen to match the s t andard cond itions in the 
BLR (IFerlandl l2006t iKorista & Goadl 120001 12004|) . Assuming 
to have solar chemical abundances and constant hydrogen den- 
sity, and using the code's AGN template for the incident con- 
tinuum shape (a continuum similar to the typical radio quiet ac- 
tive galaxies), we compute an emission line spectrum for the co- 
ordinate pair of hydrogen gas density «H[crrT 3 ] and hydrogen- 
ionizing photon flux Oh[chi V 1 ]. 

The grid dimensions span 6 orders of magnitude in each di- 
rection and, with an origin in log«n = 8, log<J>H = 17 and 
a step of 0.2 dex increments, produced a total of 961 simula- 
tions per grid. The column density A^crrr 2 ] is kept constant 
in producing the grid of simulations. Many authors claim that 
the most probable value of the column density in the BLR is 
N H _ 10 23 cirT 2 dDumont. Collin-Souffrin & Nazaroval 119981; 
IKorista & Goadl l2000l |2004|) . Since the column density de- 
fines the optical thickness of the region, here we produce 3 
grids of models changing the column density between A^h = 
[10 , 10 22 , 10 23 ]crrT 2 . Therefore, the total number of simula- 
tions is 2883. 

The output of the code includes all lineal formed in the emit- 
ting region, with a given chemical composition, hydrogen den- 
sity, and hydrogen-ionizing photon flux. We consider in our anal- 
ysis the strongest hydrogen Balmer lines (Ha to He). 

3.1.1. Modeled line ratios 

For every grid of CLOUDY models, we analyze the emission 
line fluxes by applying the BP method to the Balmer lines. In 
particular, using Eq. (3), we estimate the parameter A, from 
which we then calculate the BP temperature (hereinafter denoted 
as 7bp) of the emitting region and the error associated to the best 
fit in the BP, which we denote as /. 

A few examples of the BP applied on the simulated Balmer 
lines for a column density of A^h = 10 22 crrr 2 are presented in 
Fig. [2] In most of the cases a satisfactory fit of the Eq. (3) is 



be found at: 



- All details related to this code < 
http://www. nublado. org/ 

By default, the code gives line fluxes normalized to the H/3 flux. 
Since it has no influence in the BP analysis, we have used the normal- 
ized values. 
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Fig. 3. Contour plots for the error of the BP fit / in the ionizing flux and hydrogen density plane for different column densities 
Nh e [10 21 - 10 23 ] cirr 2 . The curves show the regions where / is smaller than 10%, 20% and 30%, while open squares represent 
the simulations for which / < 10%, and full circles are the simulations for which / < 10% and Tbp < 30, 000 K. 



not obtained and / has pretty large values. This is more notice- 
able in the plots of / in the hydrogen density and ionizing flux 
plane for all 3 grids of models, illustrated in Fig. [3] It can be 
seen in Fig. [3]that the input parameters cDh, «h for which / was 
smaller than 10% are well constrained (open squares in Fig. |3), 
and in a similar range for different column densities. The realis- 
tic models for which Tbp < 30, 000 K (denoted with full circles 
in Fig. [3) have a very high ionization parameter (low density and 
high input ionization flux), except for a couple of cases with high 
densities appearing at higher column densities. 

3.2. Theoretical line ratios 

We consid er also theoretical recombinatio n Balmer lines cal- 
culated bv iHummer & Storevl (1 19871 11992|) : IStorev & Hummer! 
(119881 U995alT We want to analyze the hydrogen Balmer lines 
obtained only by using the recombination theory, and compare 
them to the line ratios obtained by CLOUDY code, which, on the 
other hand, treats simultaneously millions of lines and physical 
processes, so that the behavior of the hydrogen lines might be 
hid den. 

I Storey & Hummerl (Il995al) calculated line emissivities, to- 
gether with other line transition parameters for hydrogenic 
atoms, assuming both Case A and Case B recombination, and 
including full collisional effects for a considerably large range 
of temperature, density, and principal quantum numbers. Their 
calculated data are given in tables av ailable in an online cata- 
logueP] (see lStorev & Hummerl 1 1995allbl for details) 

Here we use their theoretical hydrogen Balmer lines (He 
to He), calculated for the Case B recombination with the den- 
sity log « e in the range (8.2-14, 0.2 dex increment), and the 
electron temperature lying within the range 1000 to 30,000 K 
(step 1,000 K). Online data provide line emissivities within this 
range, but o nly for a small number of values of density and 
temperature (IStorev & Hummerl Il995al) . for others we interpo- 
lated line emissivit i es usin g the interpolation code provided by 
IStorev & Hummerl (Il995al) . For the BP method we use the ra- 
tios of Balmer lines, which we could directly derive from line 
emissivities. 

We plot the error of the BP fit / for theoretical Balmer line 
ratios in the electron temperature and electron density plane 
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Fig. 4. Contour plot for the error of the BP fit / in the electron 
temperature vs. electron density plane for theoretical Balmer line 
ratios. The curves show the regions where / is smaller than 10%, 
20% and 30%, while points have the same meaning as in Fig. [3] 



(Fig. 3J. Here we do not have the input ionizing flux which is 
used as an input parameter for modeled line ratios. But we can 
say that basically the electron temperature gives some indica- 
tions about the input ionizing flux. For the theoretical line ratios, 
we obtain that the / < 10% area is also covering a constrained 
area of parameter space (open squares in Fig. 0, which is resem- 
bling the one obtained from the simulated line ratios. Here the 
cases with Tbp < 30,000 K (denoted with full circles in Fig.|U) 
appear for high densities ~ 10 13 ~ 14 cirT 3 . 

The BP temperature as a function of the electron tempera- 
ture is plotted in Fig. [5] for cases when the BP fit error is less 
than 10%. Clearly there is a relation between these two temper- 
atures, especially for the electron density of 10 I4 crrT 3 (open cir- 
cles in Fig.|5]l. The BP temperature, that represents the excitation 
temperature is higher than the electron temperature. The linear 
best-fitting gives the following relation 7bp = 3133 + 1.06 x T e . 

The Balmer decrement when / < 10% follows the behavior 
of the Boltzmann-plot (solid line in Fig. |6j. For these cases in 
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Fig. 5. BP temperature plotted against the electron temperature 
for the case when / < 10%. Open circles are for the electron 
density of 10 14 cirT 3 , open squares for 10 13 < n e < 10 14 cm -3 
, and full circles for n e < 10 13 cirT 3 . The linear best-fitting of 
the cases with n e = 10 14 cirT 3 is given with the dashed line. The 
solid line marks the T BP = 30, 000 K. 
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Fig. 6. Ha/H/3 ratio plotted against the BP temperature for / < 
10%, where full circles denotes Tbp < 30,000 K and open cir- 
cles 7bp > 30, 000 K, and are separated with the vertical dashed 
line. The horizontal dashed line represents Ha/H/3 = 3.1. The 
solid line gives the line ratio calculated using Eq.(l). 

Fig.|6]the Ha/H/3 ratio is shown against the BP temperature. The 
solid line represents the Balmer decrement calculated for differ- 
ent temperatures using Eq.(l), and it is shown only for temper- 
atures smaller than 30,000 K. The averaged Ha/H/3 is 2.96 for 
cases when T^p < 30, 000 K (full circles in Fig. [6j. 

4. Observed line ratios 

We explore the usage of the BP method on the observed data. 
For that we consider measurements performed on a sample 
of 90 broad-line-emitting AGN, taken from the SDSS spec- 
tral database. For this sample the Balmer line parameters were 
accurately estimated considering the complexity of BEL (see 
lLa Mura et al.L 120071 for details). We repeat here that the line 
fluxes are corrected for the Galactic extinction, as well as for the 
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Fig. 7. Ha/H/3 ratio vs. the BP temperature for the sample of 
SDSS galaxies for objects that have the BP fit error / < 10% 
(upper panel, full circles), objects with 10% < f < 20% (middle 
panel, open squares), and with 20% < f < 100% (bottom panel, 
full triangles) are given. The dashed line represents Ha/H/3 = 
3.1. 



host galaxy contribution, which is to some extent correcting for 
the internal dust extinction, but not for the dust effects within 
the BLR. The BP method, applied to the broad component of 
Balmer lines of this sample (Fig. [TJ, produces reasonably good 
fit uncertainties, with / < 10%, in approximately 35% of cases. 
The BP temperatures are in the range T%p = 10, 300 - 25, 200 K. 
We used these cases for which / < 10% to calculate T e assuming 
that the electron density is n e = 10 14 cirT 3 and using the linear 
relation from Fig. [5] We obtain the electron temperature in the 
range T e = 7, 350-22, 250 K. The obtained values for both tem- 
peratures are in a good agreement with the values considered to 
be typical for the BLR plasma. 

We plot the Ha/Hj3 ratio against the BP temperature for this 
sample of 90 SDSS galaxies (Fig. 13 denoting differently objects 
with different /. Objects for which the BP method is working 
(full circles, upper panel Fig. |7) tend to have higher Ha/H/3 ra- 
tio (averaged value is 3.68), in contrast to objects with high fit- 
ting error / (full triangles, bottom panellT) that mostly show low 
Balmer decrements, below 3.1 (dashed line, Fig . |7J . 

One indication of the Balmer lines optical depths can be seen 
if we compare their normalized-line profiles. Some examples of 
the line-profiles comparison from the SDSS sample of galaxies 
for which / < 10% is plotted in Fig. [8] The Ha, H/3 and Hy line 
profiles (only the broad component) are compared after being 
smoothed according to their spectral S/N ratio and normalized 
to unity. From the sample of 90 objects, objects with small BP 
error / contain Ha, H/3 and Hy lines which are showing similar 
line profiles (Fig. [8), indicating that these lines may be optically 
thin. One should have in mind the relatively low resolution of 
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Fig. 8. Some examples of the line-profiles comparison for the SDSS galaxies for which the BP fit error is less than 10%. The Ha 
(solid line), H/3 (dotted line), and Hy (dashed line) broad components are shown. The name of the SDSS galaxy and the BP error / 
are given. Lines are smoothed according to the spectral S/N ratio, normalized to unity, and converted to the velocity scale. 



the SDSS, which is especially seen in the case of lines with low 
FWHM, whose line profiles are significantly affected by this. 

5. Discussion 

According to the simulation results there is a constrained param- 
eter space, at different column densities Nh = 10 21 - 10 23 cm" 2 , 
where the BP produces reasonable results with small fit uncer- 
tainties (/ < 10%). Here the cases with T BP < 30,000 K (de- 
noted with full circles, Fig. [3} that represent less than 9% of 
total number of simulations, are mostly falling within an area 
with higher ionizing fluxes and lower hydrogen densities, indi- 
cating high ionization parameters. These situations are extreme 
ones, thus the most favorable cases for the BP should happen for 
higher density, as it is seen in the case of higher column densities 
(Fig. IS. 

A similar trend is seen in case of the theoretical line profiles 
(Fig.|4]i. Here the parameter range where the BP method is work- 
ing is constrained and higher densities are preferred ~ 10 14 cm~ 3 
(Fig. E]). In these situations, even if photoionization still acts as 
the main heating source for the BLR, a contribution of the colli- 
sional processes may become relevant (especially for such high 
densities) and significantly affect the distribution of the excited 
energy configuration of the ions. The BP method reveals the oc- 
currence of such circumstances and it provides an estimate of the 
excitation temperature, which basically describes the Boltzmann 
equation of the hydrogen energy levels. 

In the sample of 90 broad-line-emitting sources from SDSS, 
there are more than 30% cases for which the BP can be ap- 
plied and it gives reasonably good fit uncertainties (/ < 10%), 
we denote these object as the "BP objects". From the observa- 
tions we see that the BP objects have higher Balmer decrement 



(Fig. with the average value of HayHj8=3.68. This value is 
higher than t he averaged value o f 3.06 of the broad-line Balmer 
decrement of iDong et al.l (120081) obtained for blue type 1 AGN. 
Our measured Balmer decrement can indicate two situations: 
(i) either there is higher contribution of collisional processes, 
i.e. we have high density plasma in the BLR (that is also sup- 
ported by the recombination theory and mod eled line ratios) 
which has been suggested in some earlier works (|yan GroningenI 
11987b ISivron & Tsurutal [T99l iBrotherton et al.l Il994l) . (ii) or 
that such high Ha/H/3 ratio is caused by the intrinsic dust ex- 
tinction effects (see discussion in §2.1). On the other hand, the 
obtained Balmer decrements and the temperature parameter A 
(see Eq.l) of the SDSS objects that follow the B P are not prefer- 
ring h igh continuum luminosities (see Fig. 1 1 in lLa Mura et all 
120071) . If we put this in the context of the modeling and theory 
findings, we can say that these objects then tend to have higher 
density in the BLR. 

We further discuss the possible mechanism of line produc- 
tion in such BLR plasma, especially the requirements to have 
optically thin Balmer lines. We start from the basic assumption 
that we have higher opacities in the Lyman line series while 
the Balmer line series is optically thin. In such a case, the Lya 
line is scattered, while the Ly/3 is almost completely suppressed. 
Therefore, we constantly have atoms excited to the energy level 
n=2. Since the energy gap between the n=2 and n=3 level is 1.9 
eV, corresponding to an electron temperature of approximately 
23,000 K, we have collisional excitation to the n=3 level, but 
also to n=4. The life-times of the levels n=3,4 are shorter than 
n=2, so we have the production of Ha and H/3 line. Even if the 
life-time of the level n=2 can be also short, this level is con- 
stantly populated by Lya scattering. Therefore, the population 
of the energy levels in the Balmer line series may be described 
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with the Boltzmann distribution, defined by the excitation tem- 
perature. 

The line profiles of the Balmer lines in the observed "BP ob- 
jects" are similar (Fig. [8), indicating that these lines are close to 
being optically thin. Therefore, we can state that the objects that 
follow the BP, the so-called "BP objects" have broad-line hydro- 
gen Balmer lines that may originate from the BLR clouds of high 
density ~ 10 14 cm~ 3 , in case the broad lines are not influenced by 
the intrinsic dust effects. Moreover, these so-called "BP objects" 
tend to have broader BEL indicating that the regions where the 
BP fit works should have higher e mitter- velocities and b e closer 
to the central engine (see Fig. 9 in lLa Mura et al. 1 120071) . 

Finally, it must be noted that the BLR plasma cannot be 
generally considered to be in the equilibrium regime for which 
the Boltzmann distribution can be used, since the density of 
emitters is probably not uniform across the region. T he den- 
sity can be expressed as (lOsterbrock & Ferlandl 120061) : N m ~ 
b,„(T, n e )g m exp(-E m „/kT), where b m {T, « e ) represents the devi- 
ation from thermodynamical equilibrium. In general, photoion- 
ization models predict that b(T, « e ) + 1, but, taking into account 
the ratio of emission lines belonging to the same series, the de- 
viation from thermodynamical equilibrium might be overtaken 
if the condition b,{T, n e )/bj(T, n e ) a 1 i s met, which should be 
possib le at least for higher energy levels (lOsterbrock & Ferlandl 
2006). 



6. Conclusions 

In this paper we present a study of the Balmer line ratios of AGN 
and their usage as a tool to investigate physical properties of the 
broad line region in AGN. The main results are based on the 
analysis of the broad hydrogen Balmer lines (Hff to He). The 
lines are obtained from the photoionization models generated by 
the spectral synthesis code CLOUDY (version C08.00), calcu- 
lated using the recombination theory for hydrogenic ions, and 
measured from the sample of observed spectra of 90 broad-line 
AGN from SDSS. The Boltzmann-plot is applied on all samples 
and the 7bp is obtained. The properties of spectra with the BP 
error / < 10% and plasma conditions for their formation are 
further explored. From our analysis, we come to these conclu- 
sions: 

i) from the CLOUDY photoionization models we obtain that 
for a limited space of physical parameter («h, ^h) the physi- 
cal processes in plasma are such, that the Balmer lines follow 
the BP. This is the case for either very high ionization param- 
eter (low density and high input ionization flux) or for high 
densities at higher column densities; 

ii) the recombination theory for hydrogenic ions predicts that 
the BP may be applicable for higher densities ~ 10 14 cm~ 3 ; 

iii) from the analysis of the sample of broad line AGN from 
the SDSS database, it follows that the BP can be applied in 
~35% of cases, and the obtained BLR temperatures are in 
the range T BP = 10, 300 - 25, 200 K. 

There is a fraction of objects that follow the BP ("BP ob- 
jects") indicating that in their BLR the physical conditions pro- 
vide that the populations of the hydrogen upper energy states 
(n > 3) follow the Saha-Boltzmann distribution. These objects 
are thus interesting and should be further investigated (e.g. what 
are their UV or X-ray spectral characteristics, or what is the con- 
nection with other AGN types) and that will be the subject of our 
future work. 



Acknowledgements. D.I. would like to thank to the Department of Physics and 
Astronomy of the University of Padova, Italy for their hospitality. Also, we thank 
the anonymous referee for useful comments. This work was supported by the 
Ministry of Education and Science of the Republic of Serbia through the project 
Astrophysical Spectroscopy of Extragalactic Objects (#176001). 

Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan 
Foundation, the Participating Institutions, the National Science Foundation, the 
U.S. Department of Energy, the National Aeronautics and Space Administration, 
the Japanese Monbukagakusho, the Max Planck Society, and the Higher 
Education Funding Council for England. The SDSS is managed by the 
Astrophysical Research Consortium (ARC) for the Participating Institutions. The 
SDSS Web Site is http://www.sdss.org/ This research has made use of NASA's 
Astrophysics Data System. 



References 

Baldwin, J. A. 1997, in Peterson B. M., Cheng E, Wilson A. S., eds, IAU Colloq. 

159, Emission Lines in Active Galaxies: New Methods and Techniques, ASP 

Conf. Sen, Vol. 113. Astron. Soc. Pac, San Francisco, p. 80 
Baldwin, J. A., Ferland, G. J., Korista, K. T., Verner, D. A. 1995, ApJL, 455, 119 
Baldwin, J. A., Ferland, G. J., Korista, K. T., Carswell, R. E, Hamann, F, 

Phillips, M. M., Verner, D., Wilkes, B. J., Williams, R. E. 1996, ApJ, 461, 

664 

Bentz, M. C, Walsh, J. L., Barth, A. J. et al. 2010, ApJ, 716, 993 
Bon, E., Popovic, L. C., Gavrilovic, N., Mura, G. La, Mediavilla, E. 2009, 
MNRAS, 400, 924 

Brotherton, M. S., Wills, B. J., Steidel, R., Sargent, W. L. W. 1994, ApJ, 423, 
131 

Canfield, R. C, Puetter, R. C. 1981, ApJ, 243, 390 
Collin-Souffrin, S., Lasota, J.-P. 1988, PASP, 100, 1041 
Collin-Souffrin, S., Dumont, S., Tully, J. 1982, A&A, 106, 362 
Crenshaw, D. M., Peterson, B. M., Wagner, R. M. 1988, AJ, 96, 1208 
Davidson, K., Netzer, H. 1979, RvMP, 51, 715 

Dumont, A.-M., Collin-Souffrin, S., Nazarova, L. 1998, A&A, 331, 11 
Dong, X.-B., Zhou, H.-Y, Wang, T.-G., Wang, J.-X., Li, C, Zhou, Y.-Y. 2005, 
ApJ, 620, 629 

Dong, X., Wang, T., Wang, J., Yuan, W., Zhou, H., Dai, H., Zhang, K. 2008, 
MNRAS, 383, 581 

Ferland, G. J. 1999, in Quasars and Cosmology, ASP Conference Series, eds. G. 

Ferland and J. Baldwin, 162, 147 
Ferland, G. J. 2003, ARA&A, 41, 517 

Ferland, G. J. 2006, Hazy, A Brief Introduction to Cloudy 06.02, University of 

Kentucky Internal Report 
Ferland, G. J., Korista, K. T., Verner, D. A., Ferguson, J. W., Kingdon, J. B., 

Vemer, E. M. 1998, PASP, 110, 761 
Ferland, G. J., Savin, D. W., eds. 2001. Spectroscopic Challenges of 

Photoionized Plasmas, ASP Conf. Ser. 247. San Francisco: Astron. Soc. Pac. 
Gaskell, C. M. 2009, New Astronomy Review, 53, 140 
Gaskell, C. M., Ferland, G. J. 1984, PASP, 96, 393 
Goodrich, R. W. 1995, ApJ, 440, 141 

Griem, H. R. 1997, Principles of plasma spectroscopy, Cambridge: Cambridge 

University Press 
van Groningen, E. 1987, A&A, 186, 103 
Halpern, J. P., Steiner, J. E. 1983, ApJ, 269, L37 
Hummer, D. G., Storey, P. J. 1987, MNRAS, 224, 801 
Hummer, D. G., Storey, P. J. 1992, MNRAS, 254, 277 

Ilic, D., Kovacevic, J., Popovic, L. C. 2009, New Astronomy Review, 53, 149 
Ilic, D., Popovic, L. C., Bon, E., Mediavilla, E., Chavushyan, V. H. 2006, 

MNRAS, 371, 1610 
Jin, C, Ward, M., Done, C, Gelbord, J. 2012, MNRAS, 420, 1825 
Konjevic, N. 1999, PhR, 316, 339 
Korista, K. T., Goad, M. R. 2000, ApJ, 536, 284 
Korista, K. T., Goad, M. R. 2004, ApJ, 606, 749 
Kovacevic, J., Popovic, L. C, Dimitrijevic, M. S. 2010, ApJS, 189, 15 
Krolik, J. H. 1999, Active Galactic Nuclei: From the Central Black Hole to the 

Galactic Environment, Princeton University Press, Princeton, New Yersey 
Kwan, J., Krolik, J. 1981, ApJ, 233, L91 
Kwan, J., Krolik, J. 1981, ApJ, 250, 478 

La Mura, G., Di Mille, F, Ciroi, S., Popovic, L. C., Rafanelli, P. 2009, ApJ, 693, 
1437 

La Mura, G., Popovic, L. C., Ciroi, S., Rafanelli, P., Ilic, D. 2007, ApJ, 671, 104 
Laor, A. 2006, ApJ, 643, 112 

Laor, A., Jannuzi, B. T., Green, R. F, Boroson, T. A. 1997, ApJ, 489, 656 
Leighly, K. M., Casebeer, D. 2007, The Central Engine of Active Galactic 
Nuclei, 373, 365 

Marziani, P., Sulentic, J. W., Zwitter, T., Dultzin-Hacyan, D., Calvani, M. 2001, 
ApJ, 558, 553 



8 



D. Ilic et al,; The analysis of the hydrogen broad Balmer line ratios 



Marziani, P., Sulentic, J. W., Negrete, C. A., Dultzin, D., Zamfir, S., Bachev, R. 

2010, MNRAS, 409, 1033 

Mathews, W. G., Blumenthal, G. R., Grandi, S. A. 1980, ApJ, 235, 971 

McGill, K. L., Woo, J.-H., Treu, T., Malkan, M. 2008, ApJ, 6 73, 703 

Negrete, C. A., Dultzin, D., Marziani, P., Sulentic, J. W. 2011. larXiv:1107,3188l 
Netzer, H. 1975, MNRAS, 171, 395 

Netzer, H., Brotherton, M. S., Wills, B. J., Han, M., Wills, D., Baldwin, J. A., 

Ferland, G. J., Browne, I. W. A. 1995, ApJ, 448, 27 
Osterbrock, D. E. 1977, ApJ, 215, 733 
Osterbrock, D. E. 1978, PNAS, 75, 540 
Osterbrock, D. E. 1981, ApJ, 249, 462 
Osterbrock, D. E. 1984, QJRAS, 25, 1 

Osterbrock, D. E., Ferland G. J. 2006, Astrophysics of Gaseous Nebulae 
and Active Galactic Nuclei (2nd ed.), University Science Books, Sausalito, 
California 

Peterson, B. M. 1994, in Reverberation Mapping of the Broad-Line Region in 

Active Galactic Nuclei, 69, 1 
Popovic, L. C. 2003, ApJ, 599, 140 
Popovic, L. C. 2006a, ApJ, 650, 1217 (an Erratum) 
Popovic, L. C. 2006b, SerAJ, 173, 1 

Popovic, L. C., Mediavilla, E. G., Kubicela, A., Jovanovic, P. 2002, A&A, 390, 
473 

Popovic, L. C., Mediavilla, E. G., Bon, E., Ilic, D. 2004, A&A, 423, 909 
Popovic, L. C., Shapovalova, A. I., Chavushyan, V. H., Ilic, D., Burenkov, A. N, 

Mercado, A., Ciroi, S., Bochkarev, N. G. 2008, PASJ, 61, 1 
Popovic, L. C., Smirnova, A. A., Kovacevic, J., Moiseev, A. V., Afanasiev, V. L. 

2009, AJ, 137, 3548 
Popovic, L. C., Shapovalova, A. I., Ilic, D., Kovacevic, A., Kollatschny, W., 

Burenkov, A. N, Chavushyan, V. H., Bochkarev, N. G., Leon-Tavares, J. 

2011, A&A, 528, 130 
Rafanelli, P. 1985, A&A, 146, 17 

Rees, M., Netzer, H., Ferland, G. J. 1989, ApJ, 347, 640 
Sivron, R., & Tsuruta, S. 1993, ApJ, 402, 420 
Snedden, S. A., & Gaskell, C. M. 2007, ApJ, 669, 126 

Shapovalova, A. I., Doroshenko, V. T., Bochkarev, N. G., et al. 2004, A&A, 422, 
925 

Shapovalova, A. I., Popovic, L. C., Collin, S., Burenkov, A. N, Chavushyan, V. 

H. et al. 2008, A&A, 486, 99 
Shapovalova, A. I., Popovic, L. C., Burenkov, A. N, Chavushyan, V. H., Ilic, D., 

Kovacevic, A., Bochkarev, N. G., Leon-Tavares, J. 2010, A&A, 509, 106 
Storey, P.J., Hummer, D.G. 1988, MNRAS, 231, 1139 
Storey, P.J., Hummer, D.G. 1995a, MNRAS, 272, 41 
Storey, P.J., Hummer, D.G. 1995b, VizieR Online Data Catalog, 6064, 
Sulentic, J. W., Marziani, P., Dultzin-Hacyan, D. 2000, ARA&A, 38, 521 
Sulentic, J. W., Marziani, P., Zamfir, S. 2009, New Astronomy Review, 53, 198 



9 



